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Abstract 

Background: Herbicide resistance in weedy plant populations can develop through different nnechanisnns such as 
gene flow of herbicide resistance transgenes fronn crop species into compatible weedy species or by natural 
evolution of herbicide resistance or tolerance following selection pressure. Results from our previous studies 
suggest that sub-lethal levels of the herbicide glyphosate can alter the pattern of gene flow between glyphosate 
resistant Canola®, Brassica napus, and glyphosate sensitive varieties of B. napus and B. rapa. The objectives of this 
study were to examine the phenological and developmental changes that occur in Brassica crop and weed species 
following sub-lethal doses of the herbicides glyphosate and glufosinate. We examined several vegetative and 
reproductive traits of potted plants under greenhouse conditions, treated with sub-lethal herbicide sprays. 

Results: Our results indicate that exposure of Brassica spp. to a sub-lethal dose of glyphosate results in altering 
flowering phenology and reproductive function. Flowering of all sensitive species was significantly delayed and 
reproductive function, specifically male fertility, was suppressed. Higher dosage levels typically contributed to an 
increase in the magnitude of phenotypic changes. 

Conclusions: These results demonstrate that Brassica spp. plants that are exposed to sub-lethal doses of glyphosate 
could be subject to very different pollination patterns and an altered pattern of gene flow that would result from 
changes in the overlap of flowering phenology between species. Implications include the potential for increased 
glyphosate resistance evolution and spread in weedy communities exposed to sub-lethal glyphosate. 
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Background 

Agricultural land represents 11% of the total surface and 
36% of the arable surface of the Earth [1] and continues 
to increase in an effort to feed a growing human popula- 
tion. As non-managed and marginal habitats are con- 
verted to agricultural use to meet this need, interactions 
between cultivated crops, associated anthropogenic se- 
lection pressures, and wild plant species increases. This 
interface represents a dynamic habitat where selection 
pressures may change quickly, creating a gradient of 
stress from lethal to survivable effects that contributes 
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to adaptation and drives the evolution of tolerance and 
resistance traits. These forces may select for increased 
weediness traits in some plant species, impacting both 
wild and cultivated environments. 

Herbicide drift is one of these selection pressures and 
occurs as a result of standard herbicide application prac- 
tices near crop fields and management targets, but can 
also occur to a greater extent when proscribed herbicide 
application methods are not followed (e.g., application in 
high wind, unregulated weed control) [2]. As a result, 
sub-lethal concentrations of herbicides impact weedy 
or native plant communities at the crop-wild interface. 
The effect of any given dose of herbicide on a plant var- 
ies greatly with species. However, field and mesocosm 
tests of sub-lethal herbicide exposure demonstrate that 
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herbicide drift can affect the plant community by redu- 
cing biomass and fecundity of both weedy and native 
plant species [3-5]. While herbicides are intended to kill 
weeds within crop fields, unintentional exposure at sub- 
lethal levels may result in the loss of species in wild and 
weedy habitats adjacent to crop fields, alter patterns of 
pollen movement between sexually compatible species, 
and change the relative contribution of different species 
to the seed bank [3,6-8]. Many factors contribute to 
the potential selective impact of sub-lethal herbicide 
exposure on weedy plant communities including: the 
genetic variation present within the community, plant 
community structure, developmental stage, inherent 
inter-specific tolerance differences, and acquired resist- 
ance via gene flow or selection [9]. Many different weedy 
species have been examined for their response to sub- 
lethal herbicide exposure and studies have shown that 
this selection pressure can be sufficient to drive the de- 
velopment of herbicide resistance. For example, expos- 
ure of weedy Lolium species to sub-lethal doses of 
ACCase herbicides has been shown to increase the level 
of resistance in progeny produced by surviving plants in 
as little as a single generation with dramatic gains in re- 
sistance in three generations both through inherited 
genes [10] and through acclimation mechanisms [11] 
such as delayed germination. While direct exposure to 
field application rates of herbicides would be expected 
to select for resistance conferred by genes of major ef- 
fect, exposure to sub-lethal levels would be expected to 
select for polygenic resistance [9]. Weedy plant popula- 
tions in field boundary habitats may be exposed to both 
strong and weak selection pressures, creating a scenario 
where resistance evolution might be optimized. 

A study system where herbicide drift selection may 
occur outside of cultivated fields is the crop Canola® 
(Bmssica napus L. [Brassicaceae]) and wild and weedy 
compatible species (see [12]) that overlap in distribution 
with Canola® cultivation. In the United States, Canola® 
production occurs primarily in the upper Midwest 
states of North Dakota, Minnesota, and Montana. Since 
their commercial release in Canada in 1995 and in the 
US in 1998, two types of transgenic Canola® have be- 
come dominant in Canola® agriculture and represent 
the vast majority of planted varieties [13]. Because of 
the overlap of compatible wild species with transgenic 
varieties, there is potential for transgene gene flow and 
hybridization between the crop and weedy species as 
well as selection for naturally evolved herbicide resist- 
ance in field boundary habitats. 

The two types of transgenic Canola® most commonly 
cultivated are varieties resistant to the herbicides 
glufosinate-ammonium (Liberty Link®), and varieties re- 
sistant to glyphosate (Roundup Ready®). Glufosinate- 
ammonium is a contact herbicide that results in the 



inhibition of glutamine synthetase, resulting in disrup- 
tions to photosynthesis and leads to plant cell death 
[14,15]. In contrast, glyphosate is a systemic herbicide 
that upon contact with plant tissues is translocated 
within the plant to growing meristems. Glyphosate in- 
hibits a key enzyme, EPSPS, in the shikimate pathway 
blocking the biosynthesis of several important amino 
acids and ultimately leads to plant death [16,17]. Because 
they each have a very different mode-of-action in target 
plants, these two herbicides are often applied in rotation 
in agricultural cropping systems. In fact, rotation of dif- 
ferent herbicides is thought to delay the natural evolu- 
tion of resistant weed populations by cycling selective 
pressures on in-field weed species [18]. 

We hypothesize that herbicide drift may affect the fit- 
ness and relative competitiveness of plants in a commu- 
nity by altering the flowering phenology of sensitive 
species without altering the phenology of resistant spe- 
cies. As a result, altered flowering phenology of sexually 
compatible feral crop and weed species may contribute 
to increased gene flow and hybridization between previ- 
ously desynchronized plants, or decrease hybridization 
between previously synchronized plants [19]. In recent 
studies, we evaluated the effect of simulated drift of the 
herbicide glyphosate at a rate of 10% of field application 
levels in constructed plant communities composed of 
transgenic and non-transgenic Brassica species [19,20]. 
Observations of plants that were treated with glyphosate 
revealed that sensitive plants appeared to have a delay in 
development resulting in a change in flowering time. 
Presumably, a sub-lethal dose of glyphosate is sufficient 
to disrupt plant development without causing mortality. 
In addition, gene flow between certain Brassica spp. var- 
ieties in these experiments was significantly increased as 
a result of glyphosate drift [20]. Based on these observa- 
tions, we conducted this study to test the hypothesis that 
sub-lethal doses of glufosinate and glyphosate change 
the flowering phenology and reproductive traits in Bras- 
sica spp. 

Methods 

Plant material and treatments 

Seven different Brassica types (hereafter, varieties) were 
used in this study. These included three crop varieties of 
Brassica napus, two wild varieties of Brassica rapa L., 
and one wild variety of Brassica nigra L. and Brassica 
juncea L each. Two of the B. napus varieties were de- 
rived from a cv. Westar genetic background representing 
a single homozygous transgenic trait in glyphosate resist- 
ant Canola® (B. napus RR), and a non-transgenic segre- 
gating variety (B. napus null) [20]. The third B. napus 
variety used was the non-transgenic B. napus cv. Spon- 
sor, which was included to determine if plant responses 
to herbicide drift can be generalized to Canola® cultivars 
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with different genetic heritage. A transgenic glufosinate 
resistant variety of B. napus was not available for these 
studies. The remaining varieties included plants grown 
from seeds of two populations of B, rapa collected from 
weedy populations in Oregon and Northern California, a 
single population of B, nigra collected from a weedy 
population in Oregon, and a single population of B, jun- 
cea (PI649101), obtained from the USDA-GRIN national 
germplasm repository. The cultivated and wild species 
used here represent a portion of a hybridization complex 
between diploid {B, rapa, B. nigra) and tetraploid (B. 
napus, B. juncea) species [12]. B, rapa and B, juncea are 
sexually compatible with B, napus but represent self- 
incompatible and self-compatible modes of fertilization 
respectively. B, nigra has not been shown to be easily hy- 
bridized with B, napus [12] but shares a genome with 
the crop species. Additionally, B, nigra is frequently 
found as a weed in the production regions of the US 
(pers obs). 

Plants were seeded in 15.24 cm (6 inches) diameter 
pots in standard potting media (Seedling Mix No. 1, 
OBC Northwest, Canby, OR) and cultivated in green- 
houses at 20-30°C temperature and 16/8 hr day/night 
light regime. Two temporal replicate experiments were 
planted 2 weeks apart (June 10, 2009 and June 24, 2009) 
with variety groups randomized and rotated in position 
on separate greenhouse benches. Replicates were exam- 
ined for a total of 100 days from the day of seeding 
encompassing the termination of flowering for the ma- 
jority of plants under greenhouse conditions. Replicates 
were examined in the same greenhouse facility and 
plants were rotated in position on the greenhouse benches 
to assure environmental uniformity. Within each temporal 
replicate, 8 individually potted biological replicates of each 
variety were examined for each treatment except for 
B, nigra and B, juncea varieties, which suffered from vari- 
able germination. In replicate one, 6 biological reps per 
treatment/control were used for B, juncea while 4 reps 
per treatment and 6 reps for control were used for 
B, nigra. In replicate two, 7 replicates were used per treat- 
ment and 9 for control for B, juncea, while B, nigra had 8 
replicates for all treatments/control. As a result, temporal 
replicate one had a total of 262 plants, while temporal rep- 
licate two had 277. 

Four herbicide stress treatments were used. Treat- 
ments involved two brand-name herbicides. Liberty® 
(glufosinate-ammonium) and Roundup® (glyphosate, iso- 
propylamine salt) applied at a simulated drift level con- 
centration of 5% (0.05) and 10% (0.10) of the field 
application rate (f.a.r.) expected near Canola® agriculture: 
(glufosinate f.a.r. = 2.48 L/Ha; 0.05 = 0.12 L/Ha, 0.10 = 
0.25 L/Ha; glyphosate f.a.r. = 2.34 L/Ha; 0.05 = 0.177 
L/Ha, 0.01 = 0.234 L/Ha). Glufosinate treatments included 
ammonium sulfate in the spray mixture (3 lbs/acre) and 



glyphosate treatments included the surfactant "Preference" 
(0.5% v/v) following suggested rates. Treatments were ap- 
plied using a track sprayer (Model RC5000-100EP, Mandel 
Scientific Company, Ltd. Guelph, Ontario, Canada). After 
herbicide applications had dried, plants were placed in 
the greenhouse and arranged in a randomized design 
to minimize spatial effects. Control plants were left un- 
sprayed. Herbicide treatments were designed to simulate 
the drift of herbicides onto escaped crop and weed popu- 
lations in adjacent non-crop habitats. As development 
times are variable between the varieties, herbicide drift 
treatments were applied 4 weeks after seeding. At this 
time, the majority of the varieties were either at the 
pre-bolting or bolting stage but no varieties had initiated 
flowering. No pollinators were released within the green- 
houses, preventing unintentional cross-pollination of var- 
ieties. Non-transgenic, self-fertile varieties {B, napus and 
B, juncea) were not restricted in the development of seed 
pods (siliques). 

Data collection 

Aboveground biomass (BIO), the total number of 
flowers (FA), the number of days to bolting (BOLT), days 
to first flower (DTF), and duration of flowering (DUR) 
were recorded for each individual plant. Days to first 
flower was recorded for all plants when the first flower- 
like structure with four petals was produced. Duration 
of flowering was recorded as the time from first flower 
to the termination of flowering (last fully formed flower) 
under greenhouse conditions. At the conclusion of flow- 
ering, plants were watered for 7 days before harvest to 
allow any developing siliques to elongate. At harvest, the 
number of flower attempts was counted by manually 
counting the siliques and pedicels on each raceme ex- 
cept for B, nigra due to the extremely large number of 
flowers on each plant of this species. Total aboveground 
biomass was collected and weighed after being dried in a 
60°C drying oven (Blue M Model POM-326E, Thermal 
Product Solutions, New Columbia, PA) for 5 days. 

Herbicide drift exposure could alter a plants ability to 
produce seeds either by impacting male function, female 
function, or both. For self-fertile species {B, napus, B, 
juncea), we evaluated the impact of herbicide treatments 
on reproduction by measuring the proportion of suc- 
cessful siliques vs. unsuccessful siliques. Measurements 
of successful self-fertility cannot distinguish reductions 
in reproductive fitness that arise either due to impacts 
on the stamen or on the pistil. Additionally, B, rapa and 
B. nigra varieties in this experiment are self-incompatible 
so additional measures of male and female function were 
conducted. Herbicide effects on male function were evalu- 
ated by digital photography and image analysis of anther 
morphology. Anthers were collected from the stamens of 
all varieties in all treatments from at least three flowers 
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per plant, and three plants per treatment. Twenty-one 
days after herbicide applications, anthers were sampled 
from freshly opened flowers and placed in a 5% sucrose 
solution and MTT viability stain [21]. We attempted to as- 
sess pollen viability with the viability stain, however, com- 
plications with pollen extraction from the deformed 
anthers obtained from glyphosate treated plants precluded 
quantitative measures of pollen viability. Instead, we quan- 
tified morphological deformities by measuring the anther 
length (L), width (W), and the W/L ratio (R) from pre- 
pared slides. Image analysis was conducted using ImageJ 
Software [22]. 

To evaluate female function, manual pollinations were 
performed between B. napus cv. RR as a paternal parent 
and B, napus cv. Null, B, napus cv. Sponsor, and B, rapa 
OR as maternal parents. Crosses were not performed on 
B, rapa CA or B, juncea due to low sample sizes of re- 
covered flowers, nor were crosses made to B, nigra due 
to high incompatibility with B, napus [12]. Pistils were 
hand pollinated at 10 days post treatment to assess the 
viability of pistils on plants in the early stages of recov- 
ery from herbicide drift. At 21 days post treatment, a 
second evaluation of pistil function on the same plants 
was conducted. The second evaluation corresponded to 
the time at which "recovered" flowers were observed. At 
least 3 individual flowers were pollinated on at least 
three plants in each treatment. Due to limited available 
pistils on B, napus plants at both pollination time points, 
it was necessary to pool the manual pollinations for cv. 
Null and cv. Sponsor varieties. The percent of successful 
manual pollinations was used to determine the viability 
of pistils at both the pre-recovery (10 day) and post- 
recovery (21 day) time points. 

Data was initially analyzed as multivariate data with 
MANOVA but due to a lack of correlation between re- 
sponse variables (data not shown), data were further an- 
alyzed with ANOVA (PROC GLM) using SAS 9.2 (SAS/ 
STAT). The two different herbicide types were examined 
using contrast statements for comparisons to control. 
Our experimental factors included Treatment (T), Var- 
iety (V), and Rep (R); all interaction effects were tested 
and included TxV, TxR, RxT, and TxVxR. When interac- 
tions were significant, examination of the simple treatment 
effects was performed [23]. Pistil viability measurements 
were analyzed using a nonparametric Mann- Whitney 
Wilcoxon Test inR [24]. 

Results 

Significant interactions between main effects were ob- 
served (Additional file 1: Table SI) indicating varieties 
should be examined separately. A significant glyphosate 
X variety interaction was expected due to inclusion of 
the glyphosate resistant B, napus cv. RR. The second 
temporal replicate had significantly longer average days 



to flower, shorter duration of flowering, reduced number 
of flowers per plant and lower biomass than temporal 
replicate one for most varieties (data not shown). How- 
ever, the differences between temporal replicates did not 
result in differences in the response of varieties to herbi- 
cide treatments but instead the magnitude of the effect 
of glyphosate treatment was greater in the second repli- 
cate (data not shown). Measurements from the two rep- 
licates were thus combined for analyses of treatment 
effects and varieties were examined for effects of treat- 
ment in contrast to control values (Table 1). 

Glufosinate treatments 

Plants that were exposed to glufosinate developed con- 
tact damage on vegetative tissues, observed as chlorotic 
and necrotic lesions, within the first few days after treat- 
ment (Figure la). After the initial plant damage, glufosi- 
nate treated plants resumed vegetative and reproductive 
growth without any further morphological indication 
of toxicity. 

Glufosinate treatment effects were primarily limited to 
the plant structure responses of aboveground biomass 
and a single effect on flower attempts. Glufosinate treat- 
ments significantly reduced the biomass produced by 
B. napus cv. Null (0.1; p = 0.004), B. rapa OR (0.1; p = 
0.0005), B. juncea (0.05; p = 0.04, 0.1; p = 0.02), and 
B. nigra (0.05; p = 0.0087, 0.1; p < 0.001) with the greatest 
reduction in biomass at the 0.10 drift level. The 
remaining three variety biomass measures were not 
significantly reduced though the data trended toward re- 
ductions at the 0.10 level (Figure 2). Glufosinate treat- 
ments did not have a consistent effect on any other 
plant response (data not shown). 

Glyphosate treatments 

Plants exposed to glyphosate demonstrated evidence of 
herbicide damage as stunting, deformation, and chlorosis 
of meristems after treatment (Figure lb). The develop- 
ment of inflorescence meristems was halted in all sensi- 
tive varieties. After a variety-specific time delay, the 
primary meristem and additional secondary meristems 
resumed development. Flowers that formed following 
treatment exposure were observed as deformed flower- 
like structures with shrunken, pale petals; these structures 
typically lacked stamens (Figure Ic). Pistil morphology 
appeared to be more resistant to glyphosate damage, 
and normal pistils were nearly always present on post- 
treatment flowers. 

In contrast to glufosinate, glyphosate treatments pro- 
duced significant changes in all plant responses measured. 
Glyphosate treatments reduced the biomass of the weedy 
B, nigra species at the 0.10 concentration. Glyphosate 
treatments also resulted in significantly greater flower at- 
tempts on both sensitive B, napus cultivars and at both 



Table 1 ANOVA results for plant measurements in response to glyphosate treatments separated for effects of 0.05 and 0.1 levels of glyphosate 
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Response variables: vegetative biomass (BIO) flower attempts (FA), changes in bolting (BOLT), days to flower (DTF), duration of flowering (DUR), male reproductive measures, pistil function, and self-fertility. Values in 
boldface type indicate significance at P < 0.05. Lack of a value indicates no significance and na indicates no measurement taken. 
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0.05 and 0.10 treatment levels (Table 1, Figure 2), possibly 
indicating a stimulatory effect of low levels of glyphosate 
on flower production in B. napus. Increased flower num- 
bers were not observed for other varieties. 

Effects of glyphosate were assessed on the flowering 
phenology pattern of Brassica spp. by examining the days 
to bolting (BOLT), days to first flower (DTF), and the dur- 
ation of flowering (DUR) (Figure 3). Glyphosate treatments 
significantly impacted all three of these measurements 
though not for every variety. Glyphosate treatments signifi- 
cantly delayed the days to bolting for B, rapa CA, B, juncea 
and B, nigra varieties (Table 1, Figure 3). Glyphosate treat- 
ments signiflcantly delayed flowering in all of the varieties 
except for the glyphosate resistant transgenic B, napus cv. 
RR (Table 1). Flowering delays were different for each of 
the varieties with B, rapa CA having the shortest delay, 



10.70 days at 0.05 glyphosate, and B. nigra having the lon- 
gest delay, 29.46 days at 0.10 glyphosate. The delayed 
recovery in flowering was more pronounced at the higher 
drift concentration (0.10) for all six sensitive varieties 
(Figure 3, Additional file 2: Table S2). Glyphosate treat- 
ments also signiflcantly reduced the duration of flowering 
for B, rapa OR, and B, nigra, 

Male and female reproductive attributes were exam- 
ined separately to determine if glyphosate drift toxicity 
affects male and female function differently. Glyphosate 
treatments typically resulted in deformed and shortened 
anthers that appear to be unable to properly dehisce and 
release pollen (Figure 4). Anther length was significantly 
reduced in aU varieties except the glyphosate resistant B, 
napus cv. RR variety. Anther width was less sensitive to 
glyphosate effects and significantly increased for B, rapa 
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Figure 2 Effect of herbicide applications on plant biomass (BIO) and flowering attempts (FA). Left hand axis represents plant biomass, 
riglit axis represents tlie number of flowers produced. Bars and points represent mean values combined from both replicates. Treatments 
denoted on x axis; glu = glufosinate, gly = glyphosate. Error bars represent +/- one SE. 



Londo et al. BMC Plant Biology 2014, 14:70 
http://www.bionnedcentral.conn/1471-2229/14/70 



Page 7 of 10 



B.juncea 



B.rapa CA 



3.rapa OR 



B.napus 
cv. sponsor 



B.napus 
cv. null 



B.napus 
cv. RR 



13^ * 



□ Bolting 
^ ■ DTF 

□ Duration 



13^ ^ 



13^ * 



O.lOGly 
0.05 Gly 
Control 



0 10 20 30 40 50 60 70 80 90 100 1 10 
Days From Seeding 

Figure 3 Changes in Bolting, Days to Flower (DTF), and 
Duration of flowering on Brassica resulting from glyphosate 
applications. Treatments indicated for each variety, ordered as 
labeled for B. nopus cv. RR Error bars indicate +/- one standard 
error. Asterisk indicates significant change in DTF at P < 0.05. Vertical 
line at 28 days indicates time of glyphosate application. 



OR, B. juncea, and B. napus cv. RR varieties. Conse- 
quently, the anther ratio was significantly different in all 
varieties (Table 1) at the 0.10 treatment level. 

Pistil function was sensitive to glyphosate drift treat- 
ments. Pre-recovery pistils had significantly reduced func- 
tion for B, napus at the 0.10 treatment level and both 0.05 
and 0.01 treatment levels for B, rapa (Figure 5). Pistils that 
were pollinated after plants appeared to have resumed 
normal flowering had much higher function in both B. 
napus and B, rapa, though function remained lower than 
pollinations made on control plants (Figure 5). 

The ability of plants to self-fertilize was examined on 
the two sensitive B, napus cultivars and the B, juncea 
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B. napus 
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B. rapa 



0.10 



I 10 days post treat 
□ 21 days post treat 

Figure 5 Effect of glyphosate treatments on pistil function, 
evaluated by successful hand pollinations. Significant differences 
noted with an asterisk. Error bars indicate +/- one SB. 



variety. All three varieties were similar under control 
conditions, producing approximately 49% ± 1% of flowers 
as siliques. The proportion of flowers that successfully 
formed a silique was significantly lower with glyphosate 
treatment, with a reduction of approximately 50% for 
aU three varieties in both 0.05 and 0.10 treatment levels 
(Figure 6). 

Discussion 

Exposure of Brassica species to sub-lethal herbicide re- 
sults in changes in biomass, flowering phenology, and 
reproductive function. Of the two herbicides tested here, 
only glyphosate exposure resulted in changes in flower- 
ing time and reproductive function. Plants that were ex- 
posed to sub-lethal glyphosate demonstrated variable 
delay in flowering time and all delays were significantly 
greater than unsprayed or glyphosate resistant plants. 
Male reproductive function was much more sensitive to 




Sponsor 



B.juncea B. nigra 



B. rapa 



Null 

B. napus B. rapa 

Figure 4 Effect of glyphosate treatments on anther length. Significant decreases in length denoted by asterisk at P < 0.05. Digital photo of 
anthers from an untreated a) B. nopus and b) B. rapa and anthers c), d) from 0.10 glyphosate treatments. 
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cv. Null cv. Sponsor 

Figure 6 Assessment of self-fertility changes due to glyphosate 

treatments. Percent of successful siliques produced on self-fertile 

varieties. Asterisk indicates P < 0.05. Error bars represent +/- one SE. 
\ J 

glyphosate exposure than female function and as a re- 
sult, plants were rendered functionally out-crossing for a 
significant period of time (2-4 weeks). These results 
demonstrate the potential for sub-lethal glyphosate to 
alter flowering time in glyphosate sensitive plant popula- 
tions. These changes in phenology in wild and weedy 
plants could contribute to changes in gene flow patterns 
between resistant and sensitive plants such that resist- 
ance alleles may unidirectionally move from resistant 
plants into functionally out-crossing sensitive plants. 
Changes in the flowering phenology and reproductive 
strategy of plants, specifically in feral conventional crops 
or sexually compatible weeds, could have important im- 
plications for transgene confinement and management. 

The results of our study demonstrate the differential 
effects of sub-lethal herbicide exposure and highlight 
the potential for ecological and evolutionary impacts in 
weedy plant communities. Evolution of herbicide resist- 
ance in weedy plant species is perhaps the greatest con- 
cern in regards to weed management [9]. Currently 
there are 221 different species that are considered herbi- 
cide resistant and weeds have evolved resistance to 152 
different herbicides [25]. Evidence for evolution of resist- 
ance has been observed for both major gene traits as 
well as multigenic traits. Direct exposure to herbicide 
spray, such as would occur within crop fields, is ex- 
pected to favor evolution of major effect genes that rise 
to fixation in the weed population quickly [26]. In con- 
trast, low and sub-lethal exposure would act to favor re- 
sistance traits that include many different loci that could 



combine to increase resistance following outcrossing be- 
tween surviving plants [10,27]. For example, the evolu- 
tion of glyphosate resistant Lolium spp. in vineyards and 
orchards [28] following repeated, non-lethal exposures. 
In this study, we observed a potential interaction be- 
tween sub-lethal exposure and herbicide resistance in 
weedy plant species. The implications of this interaction 
are that direct exposure to glyphosate would favor feral 
resistant crop plants in weedy communities, suppressing 
growth and survival of sensitive varieties while increas- 
ing the representation of resistant seeds in the weed seed 
bank. Sub-lethal exposure may also enhance the move- 
ment of transgenic resistance traits between plants 
through the synchronization and de-synchronization of 
flowering periods between resistant and sensitive species, 
creating a window of optimal out-crossing. Selective 
sterility of male tissues, but partially/fully functional fe- 
male tissues in sensitive species would explain the re- 
sults of our previous study [20] where outcrossing rate 
was seen to significantly increase for non-transgenic 
Canola® varieties exposed to sub-lethal glyphosate treat- 
ments. Expanded upon further, this temporary enhance- 
ment of pollen based gene flow between resistant crop 
varieties and sensitive weed varieties might result in the 
increased production of hybrid seed on receptive weed 
plants, impacting the structure and identity of the future 
weed seed bank. These hybrid seeds may then germinate 
and have a selective advantage in subsequent generations 
exposed to herbicide drift, contributing to the preserva- 
tion of resistance alleles in the weed population. 

Several other studies have uncovered results similar to 
ours in regard to the transient and specific loss of func- 
tion of the male reproductive structures suggesting that 
sub-lethal glyphosate effects are not unique to Brassica 
species. Studies in morning glory have revealed popula- 
tion variability in the survival of glyphosate application 
and surviving plants often have functionally female flowers 
due to abnormal stamens [29]. Similarly, studies in cotton 
have shown glyphosate-induced changes in microtubules 
in anthers, leading to poor dehiscence [30], possibly simi- 
lar to the mechanism contributing to the reduction in an- 
ther dehiscence we observed. Glyphosate resistant corn 
and cotton varieties that have reduced transgene expres- 
sion in male tissues also suffer from shortened anther fila- 
ments leading to reduced pollen transfer between anther 
and pistil [31] and reduced pollen viability [32]. Interest- 
ingly, early studies of glyphosate s mode of action demon- 
strated the function of glyphosate sprays used as a male 
specific gametocide for preventing the self-pollination of 
wheat cultivars [33] . It appears that while the utility of gly- 
phosate application for male gametocidal action is well 
known, the implications of this effect regarding gene flow 
in the environment between compatible species remain 
understudied. Though additional studies are needed, it is 
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likely that sub-lethal glyphosate exposure has the potential 
to alter the flowering phenology and mating system func- 
tion of many different wild and weedy species. 

Future studies are necessary to evaluate and describe 
the level of herbicide drift occurring in weedy plant popu- 
lations. While data exists on rates of herbicide drift under 
prescribed best practices [2], less data are available that 
describes the rates of non-regulated herbicide exposure 
and applications under adverse conditions (e.g., windy 
conditions). Additionally, little attention is paid by the ma- 
jority of weed evolution studies on weeds that grow just 
beyond agricultural fields. Instead, it is assumed that dir- 
ect exposure to herbicides is the dominant selection pres- 
sure contributing to herbicide resistance evolution. Field 
studies including multi-species plots, exposed to varied 
herbicide levels over different developmental stages would 
further add refinement to the potential implications of 
sub-lethal herbicide exposure. 

Conclusions 

In conclusion, we argue that sub-lethal herbicide expos- 
ure outside of fields may contribute to the rise of resist- 
ant weeds and our study demonstrates the potential 
mechanism for such resistance evolution. Our results 
demonstrate that sub-lethal exposure to these two herbi- 
cides results in different potential for population level 
impacts. Namely, populations exposed to sub-lethal gly- 
phosate may experience changes in flowering phenology 
that may lead to altered rates of inter and intra-specific 
gene flow. As a result of repeated exposure, it is possible 
that resistance could evolve via selection on standing 
variation in weed populations or through direct transfer 
of transgenic resistance traits due to alterations in flow- 
ering phenology and transient male-sterility. 
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